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Abstract. We study the existence and cardinality of normal geodesies of different causal 

types on M(eisenberg)-type quaternion group equipped with the sub-Lorentzian metric. We 

1 present explicit formulas for geodesies and describe reachable sets by geodesies of different 

causal character. We compare results with the sub-Riemannian quaternion group and with 

the sub-Lorentzian Heisenberg group, showing that there are similarities and distinctions. 

We show that the geodesies on H-type quaternion groups with the sub-Lorentzian metric 

satisfy the equations describing the motion of a relativistic particle in a constant homogeneous 

^ C*| electromagnetic field. 

-I— > 

1. Introduction 

The term sub-Riemannian manifold means the triple (M, d), where M is an n-dimensional 
manifold, T~L is a smoothly varying /c-dimensional distribution inside the tangent bundle T(M) 
of the manifold M with k < n, and d is a Riemannian metric defined on T~L, i. e., a positively 
definite quadratic form. Recently the study of geometric structures, where the Riemannian 
metric d on T~L is substituted by a semi- Riemannian metric g, that is a non-degenerate indefinite 
metric, started e. g., in [21 O 01 El El [Til E2]- There is no special attribution so far for such 
kind of manifolds (M,7i,g), thus we propose to call them sub-semi-Riemannian manifolds or 
shortly ssr-manifolds. In the particular case, when the metric g has index 1, an ssr-manifold 
receives the name sub-Lorentzian manifold by the analogy to Lorentzian manifold. 

In the present article we study an example of "H-type group furnished with the sub-Lorentzian 
metric. This is an interesting example not only as an almost unique known example of sub- 
Lorentzian manifold but also because it has a precise physical meaning. In the article we reveal 
the connection between sub-Lorentzian geometry and physics of relativistic electrodynamics 
basing on the example of -ff-type quaternion group equipped with the Lorentzian metric. 
We also compare characterising features of sub-Riemannian and sub-Lorentzian geometries. 
The notion of H-type groups was introduced in [9]. It is known that Riemannian manifolds 
have applications in classical mechanics. Sub-Riemannian manifolds of step 2 (Heisenberg 
manifolds) play important role in quantum mechanics. Sub-Riemannian geodesies even localy 
behave very differently from the ones in Riemannian geometry, where the energy minimising 
motion is described by a unique geodesic. A sub-atomic particle behaves in a way similar to 
an electron which moves only along a given set of directions. There can be infinitely many 
geodesies with different length joining two points. In its turn, the sub-Lorentzian structure 
underlise the motion in an electro-magnetic field. Just like space and time emerge in special 
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relativity, the electric and magnetic fields can not be considered separately. The sub-Lorentzian 
structure absorbes both phenomena, the presence of the electro-magnetic field and the the 
space-time geometry. That is why the H-type quaternion group with a sub-Lorentzian metric 
is an interesting example to work with. 

In the introduction in order to explain the main idea we would like to mention the Heisen- 
berg group as the simplest noncommutative example of H-type (Heisenberg type) groups and 
its numerous connections with physics. The Heisenberg group is the manifold R 3 with the 
noncommutative group law 

(x,y,z) o (x',y',z') = (x + x',y + y',z + z' + ~{xy' - yx')). 

Left-invariant vector fields X = d x — \yd z , Y = d y + \xd z are obtained from the group law 
and span a 2-dimensional distribution % which is called horizontal. The horizontal distribution 
can be also defined as the kernel of the contact one- form u = dz — ^(xdy — ydx) in R 3 . The 
differential of oo is the 2-form 0, = du = —dx A dy that satisfies the Maxwell's equation 
dQ, = for the magnetic field Q, in R 3 . Let us define a Riemannian metric ds 2 R = dx 2 + dy 2 
on T~L. Then the sub- Riemannian manifold M R = (R 3 , H, ds 2 R ) is also called the Heisenberg 
group. It turned out that the geodesic equation for geodesies "f(t) satisfying the non-holonomic 
constraints £ Tl(^{t)) coincides up to a constant with the Lorentz equation of motion 
of the charged particle in the magnetic field Vt. If we change Riemannian metric on T~L on 
the Lorentzian one ds 2 L = —dx 2 + dy 2 we come to the notion of sub-Lorentzian Heisenberg 
group H^ = (R 3 , Ti, dsj^). In this case the geodesic equation for non-holonomic geodesies 
coincides with analogue of the Lorentz equation for the motion of a charged particle in the 
electromagnetic field defined by £1 and the Lorentzian metric tensor. The geodesies, metrics 
properties, and other related questions on H^ were studied in [2j O 01 [11]. The lack of the 
dimension of the horizontal distribution on the Heisenberg group does not allow to reveal the 
peculiarity of the applications in the case of the magnetic and electromagnetic fields. Therefore, 
we chose the analogue of the Heisenberg group admitting a 4-dimensional distribution, that we 
called quaternionic H-type group. This example allows also to show similarities and differences 
between the sub-Riemannian and sub-Lorentzian geometries. 

The article is organized in the following way. In Section[2]we introduce the H-type quaternion 
groups endowed with different metrics: Riemannian and Lorentzian. We also present the 
differential equations for the geodesies in both cases. Section[3]is the collection of the definitions 
related to the motion of charged particles in electro-magnetic fields. We give the explanation 
of the geometrical results from the physical point of view. Section [J] is devoted to the solution 
of geodesic equations, where we fined the explicit formulas for the horizontal and vertical parts 
of geodesies. Section [5] is dedicated to study of reachable sets by geodesies of different causal 
types and estimation of the cardinality of geodesies, connecting two different points. Section [6] 
shows a brief overview of reachable sets for H^ for the sake of comparison with obtained results 
for the quaternion H-type group. 

2. H-TYPE QUATERNION GROUPS Qr AND Ql 

We remind that quaternions form a noncommutative division algebra that extends the sys- 
tem of complex numbers. It is convenient to define any quaternion q in the algebraic form 
q = a + i\b + 12c + i%d, where i\ = i 2 = i 2 = 1112^3 = —1- The scalar a G R is called the 
real part Re q = a and a vector (6, c, d) G R 3 recieved the name pure imaginary quaternion 
and is denoted by Pug. Thus q = (Reg, Pug). With this notations it is easy to introduce 
the noncommutative multiplication * between two quaternions q± and (72 using the usual inner 
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product • and vector product x in R 3 . Namely, 

(2.1) qi * q 2 = (RegiRe^i - Pugi • Pug 2 , Regi Pug 2 + Reg 2 Pugi + Pugi x Pug 2 )- 

Notice that this structure suggests itself the analogy with a Lorentzian geometry where M 4 
consists of a time part t £ M and space part (xi,x 2 , X3) £ M 3 . The conjugate quaternion q to 
q is q = a — i\b — i 2 c — 13d. It is known that quaternion q = a + i\b + i 2 c + 13d can also be 
represented in the 4 x 4-matrix form 



al + bh + cl 2 + dl 3 
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are the basis of quaternion numbers in the representation given by real (4 x 4)-matrices. 

We introduce an H-type group whose noncommutative multiplication law makes use of 
quaternion multiplication rule. Let us take the background manifold M as M7 and define the 
noncommutative law by 



(2.3) 



(x,z) o ix'z'\ 



x' z 



+ z ' + 7; Pu (^ * x ')) 



for (x, z) and (x',z') from 



\ Here Pu(x * x') is the imaginary part of the product x * x' 



defined in (|2.ip of the conjugate quaternion x to x by another quaternion x' . The introduced 
multiplication law (|2.3p makes M = (M 7 , o) into a noncommutative Lie group with the unity 
(0,0) and the inverse element (— x, — z) to (x,z). The group law defines the left translation 
L( x , z )(x', z') = (x, z) o (x', z'). Let . . . , . . . , ^ be a standard basis of the tangent 

space T p M to M at p £ M. The basic left-invariant vector fields can be obtained by the 
action of the tangent map dL^ x ^ of L^ x>z ^ to the standard basis as dL^ x z ^-^-^ = Xi(x,z), 
dL(x,z){~&r) = Zk{x,z). Then the vector fields 
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span a 4-dimensional distribution Q, which we call horizontal. The left-invariant vector fields 
Zp = ft = 1,2,3 form a basis of the complement V to Q in the TM. At each point 



1 



ANNA KOROLKO, IRINA MARKINA 



(x,z) G 
follows 



the distribution Q(x, z) is a copy of M. 4 . The commutation relations are as 



[X ,X 1 ]——Z 1 , [X ,X 2 ] — Z 3 , [X ,X 3 } — Z 2 , 
[Xi,X 2 ] = Z 2 , [Xi,X 3 ] = —Z 3 , [X 2 ,X 3 ] = —Zi. 

Therefore, {Xq, . . . , X 3 } and their commutators span the entire tangent space TM. This 
property of the distribution Q is called bracket-generating of step 2. The Lie algebra with the 
basis {Xq, . . . , X 3 , Zi,Z 2 , Z 3 } is nilpotent of step 2. 

The horizontal distribution Q can be defined by making use of one-forms. Namely, the 
one-forms 



(2.4) 



uj 1 =dz 1 - - x ^ 1 + x 3 dx 2 - x 2 dx 3 ) = dz 1 - ^da^iix, 

oj 2 =dz 2 - \{-x 3 dx° - x 2 dx 1 + x l dx 2 + x°dx 3 



co 3 =dz 3 - ^(-x 2 dx° + x 3 dx 1 + x°dx 2 - x x dx 3 



dz 2 
dz 3 



-dx T I 2 x, 

1 T 

-dx I 3 x. 
2 



annihilate the distribution Q. Here x = (x°,...,x 3 ), dx = (dx°, . . . , dx 3 ), and dx T is the 
transposed vector to dx. Thus, Q is the common kernel of forms uj^, k = 1,2,3. Let us 
consider the external differential of the linear combination uj = Ylt=i ^k^k- We get the two- 
form that is defined in 4-dimensional space 

= A dx 1 + dx 2 A dx 3 ) + $ 2 (-dx° A dx 3 - dx 1 A dx 2 ) 

+ 3 (-<ix A dx 2 + dx 1 A dx 3 ) = ]- ^ n^dx* A dx j , 



where 
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Any vector v G Q p , p G M, is called horizontal, a vector field X tangent to Q p at each point 
p is also called horizontal. An absolutely continuous curve 7 : [0, 1] — >■ M that has its velocity 
vector j(t) tangent to <2 7 (t) almost everywhere is called horizontal curve. 



2.1. Sub-Riemannian manifold Q#. Let us define the Riemannian metric (-,-)r on t ne 
distribution Q in such a way that (Xi,Xj) = 5ij, where is the Kronecker symbol. With this 
we get the sub- Riemannian manifold Qr = (R 7 , Q, (-, with the sub- Riemannian structure 

{Q, (■,■)«)■ 

Geodesies (or normal extremals) in the sub-Riemannian geometry are defined as a projection 
of solutions of the Hamiltonian equations on the underlying manifold. Let us present the 
corresponding Hamiltonian system. We denote by T p <Qr, T*Qr the tangent and cotangent 
space at p £ Qr respectively and by TQr, T*Qr the corresponding tangent and cotangent 
bundle. Thus, if (p, A) 6 T*Qr then the restriction of A on the subspace Q(p) of T p Qr is well 
defined and making use of the inner product we define a Hamiltonian function on T*Qr by 



1 3 

H R {p,\) = -Y,{\X l ) 2 , 
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where by (■, •) we denoted the pairing between vector spaces T p Q R and T*Qr. This definition 
coincides with the definition of the norm of the linear functional A over the vector space Q{p). If 
we write p = (xq, ■ ■ ■ , z\, Zi, £3) and A = (£0, • • • , £3, Si, 02, #3) then the Hamiltonian function 
can be rewritten in the following form 

H R ( P ,x) = ^\ 2 + ^\x\ 2 \e\ 2 + e^x. 

We get the corresponding Hamiltonian system 

—gA = £ + 2^ 

-\\x\ 2 e k + \^ T I k x, fe = 1,2,3 

li 

•1 



(2.6) 
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dz 



0. 



Here \9\ = (Y^t=i ^ip 1//2 ' \ x \ = (Yli=o x i Y^ 2 - After the simplification, we obtain that B k are 
constant and 

(2.7) x = Six, i£l 4 

(2.8) z k = x T I k x, k = 1,2,3. 

The solution of these equations and detailed calculations can be found in [7j. 

2.2. Sub-Lorentzian manifold Q^. Let us change the positively definite metric (•, •)# on Q 
on the Lorentz metric (that is nondegenerate metric of index 1) (-,•)/, such that 

(2.9) (X i ,X j ) L = S ij , (X ,X )l = -1, (X i ,X i ) L = l, i = 1,2,3. 

We call the triple Ql = (K 7 , Q, (•, -)^) the sub-Lorentzian manifold or the sub-Lorentzian H- 
type group and the pair (Q, (•, is named by the sub-Lorentzian structure on M 7 . 

We define the casual character on Q^. Fix a point p £ R 7 . A horizontal vector v G Q p is 
called timelike if (v,v)l < 0, spacelike if (v,v)l > or v = 0, null if (v,v)l = and « / 0, 
nonspacelike if (t>, u)l ^ 0. A horizontal curve is called timelike if its tangent vector is timelike 
at each point. Spacelike, null and nonspacelike curves are defined similarly. The choice of 
the sub-Lorentzian metric 12.91 implies that the horizontal vector field Xq is timelike and other 
horizontal vector fields Xi, i = 1,2,3, are spacelike. We call Xq the time orientation on Ql- 
Then a nonspacelike vector v £ Q p is called future directed if (v, Xq{p))l < 0, and it is called 
past directed if (v, Xq(p))l > 0. Throughout this paper f.d. stands for "future directed", t. 
for "timelike", and nspc. for "nonspacelike". 

We would like to start the description of Ql with finding geodesies, that is by definition, 
the projections of a solution of the associated Hamiltonian system on Q^. 

We construct a Hamiltonian system with respect to sub-Lorentzian metric. Locally the 
Hamiltonian function associated with the Lorentzian metric can be defined in the following 
way: 

1 1 3 

H l {p,X) = --(X,Xq) 2 + -Y,(^X 1 



2 

i=l 



If we use the coordinates for (p, A) as in the previous subsection, then the Hamiltonian becomes 

H L (£, 0, x, z) = Vt^ + liQxfnQx + fnVLx. 
2 o 
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Here £1(9) is the matrix given by (|2.5p and r] is the matrix of the Minkowskii metric tensor. 
(2.10) 7] 



( 


-1 















1 
















1 





\ 











V 



The corresponding Hamiltonian system takes the form 
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where the paricipating matrix A = tjO, is a constant matrix of the parameters and 
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Here by symbol | ^ | is denoted the expression y/6'{ + 91 + 9j. Notice that A T = -Qr]. 

After not intricate calculations two first equations of system (|2.1ip roll up to the following 
linear system of ordinary differential equations 



(2.13) 
(2.14) 



x = Ax, 



1,2,3, 



that gives the equations for geodesies on Q^. Here the conditions (|2.14|) are derived from the 
second line of the system (|2.11[) by substituting £ from the first line of this system. The exact 
formulae for obtained geodesies see in [IT] . 



3. Electromagnetic fields 

In this section we briefly introduce the notion of an electromagnetic field in order to explain 
the relation between the motion of the charged particle in an electromagnetic field and a 
sub-Lorentzian geodesic. 

Consider the Minkowski spacetime M with the Lorentzian metric r/ in it. Let (m, a) be a 
charged particle in M of a unite charge and a constant mass m with a trajectory a. Charged 
particles create an electromagnetic field and also respond on the fields created by other par- 
ticles. An electromagnetic field in M can be described by using two 3-dimensional vectors 
and B that express electric and magnetic components respectively. Electromagnetic fileds in 
the space free of charge satisfies to four Maxwell's equations 



[ ' U) Vx^ + || = 0, Vx!-| = 0. 
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where xq stands for the time coordinate in M and permittivity and permeability are supposed 
to be constant and equal to 1. If we use the covariant formulation then Maxwell's equations 
can be written in the nice symmetric form 



(3.2) 



dF = 0, d*F = 0, 



where F is a 2-form field in 4-dimensional spacetime corresponding to antisymmetric electro- 
magnetic tensor field 



(3.3) 



a0 



( 


-E 3 


-E 2 


-E x 


\ 


E 3 





Bi 


-B 2 




E2 


-B x 





B 3 




V E 1 


B 2 


-B 3 





/ 



The operator d is the exterior derivative, a coordinate and metric independent differential 
operator, and * is the Hodge star operator that is linear transformation from the space of 2 
form into the space of two- forms defined by the metric in Minkowski space, see for instance |10j . 

While Maxwell equations describes how electrically charged particles and objects give rise 
to electric and magnetic fields, the Lorentz force law completes that picture by describing the 
force acting on a moving charged particle in the presence of electromagnetic fields, see, for 
instance [Hj. This effect is described by Lorentz equation 

^ =e „ FU = F ^ U0 , 

where U is the particle's world velocity, P = mil its world momentum, t is the proper time of 
the particle, and F is an electromagnetic tenzor field. 

Let (eo, ex, e 2 , e 3 ) be any admissible basis in M; that is orthonormal, eo responds for the 
time coordinate and (e±,e 2 ,e 3 ) for space coordinates. As it was mentioned at each point of 
M the linear transformation F can be defined in terms of the classical electric and magnetic 
3-vectors E = E\e\ + E 2 e 2 + E 3 e 3 and B = B\e\ + B 2 e 2 + -6363 at that point. Set 
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The transformation A which influences on the charged particle is often called Lorentz force. 
In order to find eigenspaces of A that are invariant subspaces of this linear transformation we 
come to the characteristic equation. 

det{A - XI) = A 4 + (|^| 2 - |^| 2 )A 2 - ■ t) 2 = 0, 



where \E*\ 2 = (Ex) 2 + (E 2 ) 2 + (E 3 ) 2 and E ■ B = ExB x + E 2 B 2 + E 3 B 3 . The algebraic 
combinations |^| 2 - |^| 2 and E ■ t are the same in all admissible frames and are called 
Lorentz invariants. If both of them are equal to zero (i.e., E and B are perpendicular and 
have the same magnitudes): \t\ 2 -\E*\ 2 = t-t = 0, then A is called null transformation, 
otherwise A is said to be regular. 

Every regular skew-symmetric with respect to the Lorentz metric linear transformation 
A: M — > M has a 2-dimensional invariant subspace V such that V P| V 1 - = {0}. There exist a 
basis, that is called the canonical basis, relative to which the matrix of regular skew-symmetric 
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linear tranformation A has the form 








5 



(3.5) 
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where 6 and e are nonnegative real values such that = fe. Now 

the eigenvalues of A are easy to calculate since the characteristic equation becomes A 4 + (5 2 — 
e 2 )A 2 — 5 2 e 2 = 0, i. e., (A 2 — e 2 )(A 2 + 5 2 ) = 0, which has the following solutions: Ai j2 = ±e and 
A3,4 = 

Definition 1. The linear transformation T: M — >• M defined by 



where A 2 = Ao A, I is the identity transformation I{x) = x for every x £ M and ti(A 2 ) is the 
trace of A 2 , is called the energy-momentum transformation associated with A. 

Observe that T is symmetric with respect to the Lorentzian inner product and is trace-free, 
i. e. trT = 0. The term — T 1 = Teo • eo = g^(|i^| 2 + |I^| 2 ) is called the energy density in 
the given frame of reference for the electromagnetic field F with the form fj3.3[) . The 3- vector 

±E x B = {E 2 B 3 - E 3 B 2 )ei + (E 3 Bi - E x B 3 )e 2 + {E\B 2 - E 2 B x )e 3 is called the Poynting 

3- vector and describes the energy flux of the field. Finally, the 3x3 matrix (1^)^=1^,3 is 
known as Maxwell stress tensor of the field in the given frame. Thus, the content of the matrix 
of T determines the energy content of the field F in the corresponding basis. 

4. SUB-LORENTZIAN GEODESICS AND THE TRAJECTORIES OF THE PARTICLES 

The noncommutative multiplication law (|2.3p of quaternion H-type group (R 7 , o) defines the 
non-integrable distribution Q = span{AQ), . . . ,^3} or in the covariant language the nonholo- 
nomic constraints uj\ = oo 2 = uj 3 = 0. The curvature of the distribution gives rise to the skew 
symmetric transformation in 4-dimensional space (|2.5p . Independently whether this space 
has Euclidean structure or it is the Minkowskii space the antisymmetric 2-form £1 defines the 
electromagnetic field, since it trivially satisfies the Maxwell equations (|3.2|) . We emphasize that 
the geometry of the nonholonomic manifold (R , o) is related to the geometry of 4-dimensional 
space where a constant electromagnetic field acts. Given a positively definite metric (-,-)r 
(Riemannian metric) on the nonholonomic distribution Q we obtain a sub-Riemannian mani- 
fold Qr. The Hamiltonian function Hr in this case is reduced to the Lorentz equation in the 
Euclidean space R given by (|2.7p . Equation (|2.7p describes a motion of charged particle of 
unit charge in magnetic field in 4-dimentional Euclidean space. Replacement of a positively 
definite metric (-,-)r on a nondegenerate indefinite metric, the Lorentzian metric (•, leads 
to the relativistic Lorentz equation (|2.13p in the Minkowskii space. It has more connections to 
physics since it is related to the motion of the charged particle in electromagnetic field in the 
Minkowskii space that is closely connected to the general relativity. 

Our aim is to find geodesies that are projections on 1Z 7 of the solutions to the Hamiltonian 
system for Hi defined on Q. The Hamiltonian system is reduced to the equations (|2.13p 
and (|2.14p . Equations (|2.13p are ordinary differential equations in 4-dimensional space with A 
that is skew-symmetric with respect to the Lorentzian metric n. This makes us to endow the 

4- dimensional space with the Lorentzian metric, producing the Minkowskii space M = (R , 77). 
Therefore, at each point of M we can elect a model of electromagnetic field which corresponds 
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to a linear skew-symmetric transformation A : M — > M that assigns to the world velocity x = U 
of a charged particle passing through that point the change in world momentum ^j? that the 
particle should expect due to the field. Since we can assume that the charge and the mass 
of the particle equals 1 we get P = U and the Lorentz force law becomes =£■ = AU that 
is exactly equations (|2.13p . If we set E\ = —B\ = 02, E% = —E>2 = 03, E% = —B3 = 0\, 
then we conclude that equation f)2 . 13|) describes the motion of a particle of unit charge in the 
constant electromagnetic field with E = — B. In this case one of the Lorentz invariant is zero: 
\B\ 2 - 0\ 2 = 0, and the other E ■ B is equal to -\0\ 2 , where \0\ 2 = 2 + 2 + Q\. The matrix 
A 2 given by (|2.12[) has zero trace. Therefore, the energy momentum transformation is given 
up to a constant multiplied by symmetric matrix A 2 . The energy density is equal to — 2\0\ 2 . 
The Poynting 3- vector in our case is zero vector. The Maxwell stress tensor is given by 

— 20x03 —10\02 
(4.1) i -2010 3 -02-02 + 02 2 2 3 




20 2 3 -0{ + 



2^ 2 2 -0l 



Let us find a canonical basis for the matrix A. First we find a plane in M 3 in which both 
of the vectors E and B are lying, then we rotate it in such a way that it will coincide with a 
plane X3 = in M 3 . That is, let us choose a right-handed orthonormal basis {ei, (22,63} of the 
space M 3 = span(ei, e 2 , 63) in which £3 = B3 = 0. There are infinitely many planes containing 
both vectors E and -E*. Fix one of these planes: for example, the one passing through the 
axis X2 = (0, 1, 0): 

(4.2) - Xl E 3 + x 3 E x = 0. 

With the help of the rotation 

(cos a sin a 
1 
— sin a cos a 

we turn it on the angle a = arccos , El between the plane (|4.2|) and x% = around the axis 

V E l + E 3 

X2 so that the third coordinate of E equals to zero. Consider now a Lorentz transformation 

Rl = ( R 

of a basis {eo, ei, e%, e^}, where R is given by (|4.3p . It yields a new admissible coordinate 

system {e , ei, e 2 , e 3 } in which E = (y/ ' E\ + E$,E 2 ,0) and B = {-\JE\ + Ej, -E 2 ,0) and 
matrix A can be defined in the following way: 
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where E x = y/E 2 + E 2 = -B x = y/0 2 + 2 , E 2 = E 2 = -B 2 = 3 , and % = -B 3 = 0. 

Since A : M — > M is regular skew-symmetric transformation, then it has a 2-dimensional 
invariant subspace U such that U Pi U = {0}. Then U 1 - is also a 2-dimensional invariant 
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subspace for A and there exist real numbers e ^ and 5^0 such that 

A 2 u = e 2 u for all u = (ux,u 2 ) G U and 
A 2 v = -5 2 v for all v = (vi,v 2 ) & . 



Take any future directed unit timelike vector e~o in U, for example (1, 0, 0, 0) = V1 ^tg! 2 , where 
t;, are eigenvectors of matrix A (see [11] ) . Then spacelike unit vector e"i in U and the real value 
e ^ can be found from the conditions Ae~o = ee~\ and Ae\ = ee*o. We get e"i = (0, — , j^, ) 

and e = |< 



C2 



V3+V4 



and Ae*3 = — 5e 2 



Now, let e2 be an arbitrary unit spacelike vector in U , for example, select 
J (0,6x0s, 9 2 + 9 2 ,— 6*26*3). Then construct £3 satisfying Ae 2 = 5e~3 



\e\VW- 



We obtain £3 = (0, 



72P3; 



■ 0, 



is an orthnormal basis for M, which is called canonical and 



and 5 = \9\. Thus, {eo, e%, e 2 , e^} 



( \e\ 



(4.4) 



A 



\ 











-\e 




/ 



Electric and magnetic fields corresponding to this transformation are E = ee*3 = \6\e3 and 
B = 5e~3 = \O\e~3, so that an observer in this frame will measure them in the same direction (of 
X3-axis) and of magnitude Matrix A is a block-type matrix, where (2 x 2)-block in the left 
upper corner coincides with a matrix in sub-Lorentzian Heisenberg case and the right lower 
(2 x 2)-block coincides with a matrix in usual sub-Riemannian Heisenberg case (see |11|). 

Let us denote by P the matrix which columns are orthnormal basis vectors eo, . . . , £3. Then 
in new basis A = P~ l AP and the vector x G M is of the form x = Px id, where x id is a vector 
in old coordinates. Therefore, x id = P~ 1 x. 

It is clear that real eigenvalues of A are A = ±\9\ from the canonical form ()4.4p . Therefore, 
eigenspaces are span{e"3 + e^} and span{e*3 — e^} respectively. Directions £3 ± £4 are called 
principal null directions of A. 

The Lorentz equation (|2. 13[) in canonical coordinates takes the form 



/ # \ 

at 
dxi 
dt 
d±2 
dt 

\ dt / 



( 

1*1 



V 













-\e 

\9\ 



\ 



/ x \ 

x 2 

V ±3 J 



( Wx \ 

1^0 

1*2 / 



V 



It splits into 2 independent systems 




Fl, 
'|*o 



^ = "1^3, 



dt 



\X 2 



We wish to solve this system under the initial conditions x(0) = 0, x(Q) 
for x(t) is 

± (t) = ±o(0) coeh(|0|*) + ±i(0) sinh(|0|i), 
= i (0) sinh(|0|i) + ii(0) cosh(|0|t), 
± 2 (t) = x 2 (0) cos(\6\t) - £3(0) sin(|0|t), 
x 3 (t) = i 2 (0) sin(|0|t) + i 3 (0) cos(|0|t). 



x . 



The solution 



(4.5) 
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Integrating these expressions, we get the solution x(t), that we write in the matrix form as 
x(t) = Wxq with 



/ sinh(|0|i) cosh(|0|t) - 1 



(4-6) W = ^ 



\ 



\ 





cosh(|0|i) - 1 sinh(|0|t) 

sin(|(9|t) cos(|0|i) - 1 

l-cos(|0|t) sin(|0|t). / 

Now, let us introduce following notation: v\ = — £q(0) + if (0), v 2 = £f(0) + £§(0). 

Lemma 1. The projection of the geodesic onto (xq, x\) -plane is a brunch of the hyperbola with 
the canonical equation 



(4.7) 

Proof. Since 



xo + 



si(0) 



X\ + 



£ (0) 



x {t) = 777(^0(0) sinh(|0|t) + £i(O)(cosh(|0|t) - 1)) 



Mt) = (x (0) (cosh 1) +xi(O)sinh(|0|t)), 

we calculate that — x^(t) + x\{t) = 2 ^ x o(o)+ a i(°)) . g^^ 2 MJ*_ This expression can be rewritten 



□ 



as stated in (j4.7j) 

Lemma 2. The projection of the geodesic onto (x2,x^)-plane is a circle with the center at 

W °f the radius 1C" 
Proof. Since 

X2(t) = 7757(^2(0) sin \9\t + £ 3 (0)(cos \0\t - 1)) 



^(t) = 7^7 (x 2 (0)(l - cos \6\t) + £ 3 (0) sin(|0|t)) , 



we get x\(t) + x\{t) = ftp • sinh 2 The expression leads to 



(4.8) 



x 2 + 



£3(0) 
1*1 



+ ^3 



£ 2 (0) 



V2 



The horizontality conditions (|2,8p in canonical basis have the form 

i fc = {IkX id,x id) = {hP~ X x,P~ x x) = (PI k P~ l x,x) 

where the matrices 



Ji = P/1P- 1 



1 



0(0), 



J 2 = PJ 2 P~ 



( 2 0i \ 

-02 01 

-0i 2 

V -01 -0 2 j 



□ 
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J3 = PUP- 1 



\o\y/w+n 



V 



#1#3 

n + el 

— 02^3 



C1C3 


#2 #3 



ol + el) 

@2@3 


6163 



P2C 3 
?? + fl| 
— 6>1#3 







/ 



are skew-symmetric with respect to usual euclidean metric. More explicitely, z = Vz, where z 
are coordinates in the canonical basis, z = (Si>S 2 , S 3 ) are auxiliary expressions z^ = (I^x, x), 
x is given by ( |4.5p and 



/ 



(4.9) 



V 



ex 

02 



02 



0103 



e 2 e 3 



83 

~¥\ 




\ 



J 



Notice that V is an orthogonal transformation in M 3 while the matrix P represents the orthog- 
onal transformation M. . It is more convenient for us to work with the expressions z\, z 2 , Z3. 
Then zi, Z2, Z3 can be obtained by the orthogonal transformation V of zi, Z2, Z3. Taking into 
account that x = Wx° and (I4.5|) . we calculate 

= j^(«l(cosh(|0|i) - 1) + «2(cos(|0|t) - 1)), 



1 

1*1 L 



2 cosh |0|icos |0|tx o (O)x 2 (O) + 2 sinh \6\t cos |6>|fcci(0)x 2 (0) 

- 2sinh|6»|tsin|(9|txi(0)i;3(0) - 2 cosh \6\t sin |6>|ix (0)x 3 (0) 

+ cosh(|0|i)(xi(O)x 3 (O) -x (0)x 2 (0)) +sinh(|e|t)(io(0)i 3 (0)-ii(0)a? 2 (0)) 

- cos(|0|i)(x o (O)x 2 (O) + ii(0)i 3 (0)) + sm(|0|i)(x o (O)x 3 (O) - ii(0)i 2 (0)) 



1 

1^1 



- 2 cosh I (9 1 1 cos I (9 1 tx (0)x 3 (0) - 2 sinh 1 <9 1 1 cos 1 <9 | t±i (0)x 3 (0) 

-2sinh|0|tsm|6»|txi(O)x 2 (O) - 2 cosh \6\t sin |0|ix o (O)x 2 (O) 

+ cosh(|0|t)(x o (O)x 3 (O) + xi(0)x 2 (0)) + sinh(|6>|i)(x (0)x 2 (0) + xi(0)x 3 (0)) 

+ cos(|0|i)(x o (O)x 3 (O) - ±i(0)x 2 (0)) + sm(\9\t) (x (0)x 2 (0) + Xi(0)x 3 (0)) 

Let us use the following notation for the constants 

01 = x (0)x 3 (0) + xi(0)x 2 (0), a 3 = x (0)x 2 (0) + xi(0)x 3 (0), 
a 2 = io(0)x 3 (0) - xi(0)x 2 (0), a 4 = x (0)x 2 (0) - xi(0)x 3 (0). 



Thus, we integrate 



h(t) = (yi(smh(\e\t) - \9\t) + u 2 (sin(|0|t) - \9\t 



(4.10) 
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z 2 (t) = |^(ai(cos(|0|t)cosh(|0|*) - l) 

-a 2 (sm(\e\t)smh(\e\t + coa(\e\t)-cosh(\e\t))) 
+ a 3 (cos(|0|i)sinh(|0|t)-sin(|0|i)) 

+ a 4 (sin(|0|i) cosh(|0|i) - sinh(|0|i)) 
z 3 (t) = ^(a 4 (cos(\9\t)cosh{\9\t) - l) 



a 3 (sin(|0|*) sinh(|0|t + cos(|0|t) -cosh(|0|t))) 
a 2 (cos(|0|t)sinh(|0|t) -sin(|0|t)) 

ai(sin(|0|t)cosh(|0|t)-sinh(|0|t)) 



Observe that 

al + a\ = al + al = (x§(0) + ±l(0))v 2 , 
a\a 2 + 0304 = —V1V2, aia 3 - a 2 a 4 = 2xo(0)xi(0)w 2 . 

Then the direct calculations yield 

n ~4(t) + z 2 3 (t) = |j!(cosh(|0|t) - cos(|0|i) - sin(|0|t) sinh(|0|t)) 

((so(0) +x?(O))cosh(|0|i) + 2x (0)xi(0) sinh(|0|i) - vi). 

5. Reachable sets by geodesics 

We wish to describe the set of points in Ql that can be reached from the origin (x, z) = (0, 0) 
by a geodesic: timelike, lightlike or spacelike. We can fix the starting point at the origin O = 
(0, 0), since the solutions of the Hamiltonian equations are invariant under the left translation. 
We start from the simple lemma, that is related to the case of null transformation A. 

Lemma 3. // \9\ = 0, then the system (|2.1ip with initial data x(0) = z(0) = 0, £(0) = £°, 
r/(0) = r/o has the solution x(t) = z(t) = 0, £(i) = £°, rj(t) = rp . The projections to 
(x, z)-space are straight lines that are timelike if \£, \ L < 0, lightlike if \C°\ 2 L = 0, and spacelike 
if £° 7. - <>. 

Proof. The condition \9\ = immediately implies that £(i) and r](t) are constant and x(t) = £°t. 
Then Ef. = §(£°) r -?fc£ = since the matrix 1^ is skew symmetric. □ 

The description of the reachable set by causal curves is very complicated in general, therefore, 
we present here some particular cases. We mostly reduce our considerations to the sets that 
can be reached by geodesics, since the geodesics can not change its causal character. 

5.1. Connectivity by geodesics between (0,0) and (x ,z l ), where \x l \\ = 0. We need 
to solve the equations 

x — A.x ^ Zfc — x J ]%x ; k — J . 2.. 3 
with the boundary conditions 

(5.1) x (0) = z(0) = 0, x(l) = x\ z{l) = z l , where \x l \\ = (r]x l ,x l ) = 0. 
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We find the relation between the initial velocity and the value x(l). Substituting t = 1 in 
x(t) = W±°, we obtain 



(5.2) 



It 1 ! 2 
F \L 



4 
W 



±1(0) + ±1(0)) sinh 2 M + (x 2 (0) + iftO)) sin 2 M 



We also have ±(0) = x° = W 1 (t)x(t), where 

/ sinh(|fl|t) 

101 

(5.3) 



w-\t) = O 

v ; 2 



cosh(|0|t)-l 

-1 



V 








-1 

sinh(|fl|f) 
cosh(|0|t)-l 









sin(|0|t) 
l-cos(|0|t) 

1 






-1 

sin(|g|t) . 
-cos(|6»|t) / 



from (|4.6p . Putting i = 1, we calculate 



|~0|2 

F II 



(rj±°,±°" 



!) + *?(!)) 



1 



sinh 10 



+ (*!(!) +*!(!)) 



sin 



The expression (|5.2p can be written in the form 







It- 1 ! 2 

F II 



\0\ 



(vi sinh z ^ + v 2 sin — ). 



(5.4) - „ |<>|2V -.— 2 ■ — 2 

Let us describe all possible initial velocities which lead to vanishing of the norm (x 1 ) 2 . 

Case 1. \0\ = 0. In this case Lemma [3] implies \x l \ 2 L = |x°| 2 and we conclude that if (x , z ) 
is such that \x | 2 = and z 1 = then there is a unique geodesic, that is lightlike straight line, 
connecting (x 1 ,^ 1 ) with (0,0). 

Case 2. v 2 = 0, \0\ ^ 0. Therefore, v\ = 0. In this case \±®\\ = and the possible 
geodesic is lightlike that remains lightlike for all t 6 [0, ±oo]. Let us write the equations 
for these geodesies taking into account the condition on the initial velocity ±o(0) = =txi(0), 
x 2 (0) = i 3 (0) = 0. 
(5.5) 



If xq(0) = xi(0), then xo(t) = x±(t) 



±o(0) 



(sinh \0\t + cosh \0\t - 1), x 2 (t) = x 3 (t) = 0, 



(5.6) 
if x (0) 



-xi(0), then xo(t) = —x\(t) 



ao(0) 
1*1 



(sinh \0\t - cosh \0\t + 1), x 2 (t) = x 3 (t) = 0. 



This is not obvious parametrization of straight lines x(s) = (xo(0)s, xi(0)s, 0, 0), Xq(0) = ±1(0). 
In both cases Zk(t) = ± T (t)Ikx(t) = 0, k = 1,2,3 that implies Zk(t) = Zk(t) = 0. We conclude 
that the origin can be joined with a point (xo(l), xi(l), 0, 0), where Xq(1) = x\(l) by lightlike 
geodesic that is straight line. Thus, case 2 is a particular situation of the case 1. 

Notice that the equations of the horizontal part (|5.5p and (|5.6p of geodesies in coincides 
up to reparametrization with the equations of horizontal part for geodesies in provided the 
initial velocity ±o(0) = ±xi(0) (see [TT]). 

The results of cases 1 and 2 may be united in the following statement 



Theorem 1. Let A = (x 1 ,^ 1 ) be a point such that \x l \\ = 0, Zk = 0, k = 1,2, 3 and 
Xq(1) + x 2 (1) 0. Then there is a unique lightlike geodesic joining the origin and A which is 
a straight line. 
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Case 3. v 2 7^ 0, \Q\ ^ 0. We write the expression (j5.4j) in the equivalent form 



(5.7) 



sin h£ v\ 



sinh 2 M " v 2 

and consider the following subcases. 

3.1) Let v\ = then |x°| 2 = v 2 > and in this case all geodesies are spacelike. Then (|5.7p 
implies that |0| = 27m, n G N \ {0}. 

3.1.1) Let us suppose that xo(0) = ±i(0) = and x 2 (0), £3(0) are arbitrary. The main result 
is expressed in the following 

Theorem 2. Given a point A = (O,^ 1 ) there are uncountably many spacelike geodesies con- 
necting the origin with A. The geodesies are given by equations 

2\z 1 \ 

x (t) = xi(t) = 0, xl(t) + xl(t) = - — - sin 2 (7rnt)) 

irn 

'sm(2nnt) 



^) = kl(=^-«)(«(D,«(i),».(i)). 



The geodesies have the lengths y / 2irn\z 1 \, n G N \ {0}. 

Proof. We start from the equations for x-coordinates of geodesies. Since x(t) = Wx° we get 



(5. 



x (t) ( \ 

xi{t) _ J_ 

X2(t) ~ 2nn x 2 (0) sin(2vrnf) + x 3 (0)(cos(2vrnt) - 1) 
x 3 (i) \ i 2 (0)(l -cos(27rnt)) +x 3 (0)sin(27rnt) J 

by (|4.6p . The projection into (x2, X3)-plane is 

(5.9) xl(t) + z§(i) = 4v 2 sin 2 (vrnf)), 

which by lemma ([2]) can be rewritten as the canonical equation of the circle on the (z 2 , Z3)-plane 
with the center at (-^, ^) of the radius ag[: 

(5.10) x 2 + — + z 3 1 - 



2-Kn J \ 2irn J (27rn) 2 

The number n G N \ {0} reflects the number of turns along the circle for any fixed initial 
velocity V2 and the radius of the circle. To exclude V2 from the equations we use ([4.10p and 
find 

h(t) = —f^(sm(2Trnt) - 2imt), z 2 {t) = 0, z 3 {t) = 0. 
47r z n z 

Set t = 1 in the expression for zi(t) and get v 2 = — 27rnzi(l). Since z(i) = Vz(t) with "P given 
by P~9|) we get 

(5.11) ,i = ?1 h)(Zgl 02 g^g— |0|=2vrn. 

It implies that —23.(1) = \z l \ and ^2 = 2im\z 1 \. Fixing n G N \ {0}, we fix the speed v 2 of a 
geodesic, but we still have the choice in the directions of (£2(0), £3(0)) that are parametrized 
by the unit circle. It gives uncountably many geodesies. 

We get the equations for z-coordinates from z(t) = Vz(t) taking into account that the values 
of Ok related to the values of z 1 by (|5.1ip . 
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Since the geodesies are spacelike the length £(7) of a geodesic 7 can be calculated from the 
formula 

-l 



1(7) 



v 2 dt = \J 2im\z 







□ 



Remark 1. Notice that in Lorentzian Heisenberg group Ml there is no geodesic of any causal 
type joining (0, 0, 0) and (0, 0, z) with z 7^ (see [TT]). 

3.1.2) We suppose now that xq(0) = xi(0) 7^ 0. We present first the auxiliary calculations 
and then formulate the main statement. We have 

ai = a 3 = x (0)(x 2 (0) + x 3 (0)), -a 2 = a 4 = x (0)(x 2 (0) - x 3 (0)). 

Then the equations of geodesies take the form 



(5.12) 



x (t) 

XI (t) 
x 2 (t) 

x 3 (t) 



( 



2irn 



x (0)(e 27rn * - 1) 
x (0)(e 27rn * - 1) 
x 2 (0) sm(2vrni) + x 3 (0)(cos(2vrnt) - 1) 
\ x 2 (0)(l - cos(2vrnt)) + x 3 (0) sin(2vrnt) J 



The relations P~TU]) lead to 



(5.13) 



v 2 



z 2 (t) 



*o(0) 
(2vm) 2 



x 2 (0)((cos(2™i) - 1) (e 2nnt + 1) + sin(2vrnt)(e 27rn * - 1] 
+x|]((cos(27rnt) + l) (e 2 ™* - l) - sin(2vrnt) (e 2nnt + 1) 



i (0) 



(27m) 2 



£ 2 (0) ( ( cos(27rni) + 1) (e 2 ™* - l) - sin(2vrnt) (e 27m * + 1 
-x 3 (0)((cos(2vrnt) - l) (e 2lTnt + l) - sin(2vrnt) (e 2wnt - 1 
where n G N \ {0}. Setting t = 1 in (|5.12j) we get 



(5.14) 

The equations (|5.13p imply 
(5.15) v 2 = -5i(l)2vrn, 



x (0) =x (l) 



2irn 



2nn 



x (0)x 2 (0) 
(2vrn) 2 



2 e 



27rn 



1) 



x (0)x 3 (0) 
(2vm) 2 



2(e 



27rn 



1) 



Theorem 3. Given a point A = (x 1 ,^ 1 ), \x 1 \l = 0, xo(l) = xi(l) 7^ ; \z l \ 7^ 0, there are 
uncountably many spacelike geodesies connecting the origin O = (0, 0) with A. The equations 
of geodesies are given by expressions 



(5.16) 



x (t) = xi(t) = x (l)- 



,2nn 



x 2 {t)+x 2 (t) 



2*1 (1) . 2 



TTTl 



sin (irnt)) 
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and z(t) = Vz{t), where 

Zl (t) = Zl (l)(t 

~Z2it) = ^ _ [S 3 (l)((cos(27rnt) - l)(e 2 ™< + 1) + S in(2vrnt)(e 2 ^ - 1)) 

+5 2 (l)((cos(27mi) + l)(e 27rn * - 1) - sin(2™i)(e 27rn * + 1))] , 
W = 2{e2 l _ 1} [%(l)((cos(27rni) + l)(e 2 ™< - 1) - S in(2vrnt)(e 2 ^ + 1)) 

-5 2 (l)((cos(27rnt) - l)(e 2 ™* + 1) + sm{2Trnt){e 2nnt - 1))] , 
V~ X z x , and n € N\ {0}. 



Proof. Fix n 6 N and choose any 81,62, #3 such that y 5Zfe=i l^lfe = 27rn. We have uncountably 
many triples that are parameterized by the 3-d sphere. This choice defines the orthogonal 
transformation V given by (|4.9p . Given z l we find the auxiliarry parameters z 1 by z 1 = V z . 

Setting the values of io(0), t> 2 and x 2 (0),X3(0) from (|5.14p and (|5. 15[) into the general 
solutions (|5.12|) and (|5.13|) and applying the orthogonal transformation V to z(t) we finish the 
proof. □ 

We observe that projections of the geodesies onto (xq, xi)-plane are straight lines xq(s) = 
xi(s) = s. The projections onto (x 2 , X3)-plane are circles from the lemma [2j 

3.1.3) For the case xo(0) = — xi(0) 7^ in the same way as above we obtain the following 

Theorem 4. Given a point A = (x 1 ,^ 1 ), \x 1 \i = 0, Xq(1) = — X\(l) 7^ 0, \z 1 \ 7^ 0, there are 
uncountably many spacelike geodesies connecting the origin O = (0,0) with A. The equations 
of geodesies are given by expressions 

g—2nnt j 

x (t) = -xi (i) = x (l) • 



and z(t) = Vz(t), where 



= — 27T71 I 

x\(t) + xl(t) = ~ 2zi ^ sin 2 (7rnt) 
Tin 



~ / s ~ /-w sin(27rnt) , 
Zl(t) = z 1 (l)(t- 



2im 

m = 2(e - 2 l + 1) [-%(l)((cos(2vrnt) - lXe" 2 ™' - 1) - sin^TmtXe" 2 ™' + 1)) 

+5 2 (l)((cos(2vrnt) + l){e~ 27Tnt + 1) + sin(27mt) (e~ 27mt - 1))] , 

hit) = 2(e _ 2 ! - [5 3 (l)((cos(27rnt) + l)(e'^ nt + 1) + sm^Ke" 2 ™' - 1)) 

-5 2 (1)((1 - cos(2Trnt))(e- 2 ™* - 1) + sin(27rni)( e - 2 ™' + 1))] . 

z 1 = V~ X z x , and n G N\ {0}. 

3.2) Consider the case v\ > 0, i. e., the initial velocity x° is spacelike. It suggests that the 
ratio ——is negative, which contradicts to (|5,7p . Hence, in this case there are no spacelike 
geodesies with the property (|5.ip . 



3.3) Now suppose v\ < 0. It implies that the ratio — ^ is positive and here we have 3 more 
subcases: 
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3.3.1) — ^ > 1. In this case the initial velocity appears to be timelike. There are no points 

■ 2 l» 

sin 1 — L 



where this can happen because the function v(\9\) = ^ 2 ^ is wave decaying with exponential 

speed on the right halfplane and bounded with from below and with 1 from above. Thus, 
there are no timelike geodesies with the considered timelike initial velocity. 

3.3.2) — ^ = 1. Then sin 2 y = sinh 2 -y-, which is possible if and only if |0| = 0. This case 
corresponds to Case 1), when the initial velocity is lightlike. 

3.3.3) < — ^ < 1. Let us investigate the quantity of Hamiltonian geodesies in this case. 
Using (pETT) and the fact that x° = W~ 1 (1) x 1 we calculate that 

(5.17) i 2 2 (l) = T^Pfw ( C ° sh l g ' " COS l g ' " Sin l g ' Sinh l g D ' 

2\9\ z sinh ^ sin ^ 

= ^^T^TW ( sill2 T( SillIl l g l - 1^1) "Sinh 2 ^(sin |^| - \6\)) 
4|0| 2 sinlr ^ sin 2 1 1 

and notice that the ratio ^-Hl^i/ll depends only on 101 and 

-2 Zl (1)^(1) ^ 

z 2 (l) + z 2 (l) cosh 101 — cos 101 — sin 1 1 sinh 1 1 

(5.18) 



2Z 1 (1)X 2 (1) sin 2 M . ( sinh \Q\ _ _ sinh 2 M . ( sin \Q\ 

functior. 
(5.19) n(\9\) 



2 

Lemma 4. T/ie function 

cosh |0| — cos |0| — sin \9\ sinh |0 



sin 2 M • (sinh |0| - |0|) - sinh 2 M • (sin |0| - |0|) 

is nonnegative, has countably many points \9\ k where fi(\9\k) vanishes and on each interval 
{\9\k, \9\k+i) > k = 0, 1, 2, . . . the function fj, has a unique critical point m^. On each of these 
intervals the function [i strictly increases from to //(m^), and then, strictly decreases from 
yu(mfc) to 0. 

Proof. Let us denote the numerator cosh \ 9\ — cos \ 9\ — sin \ 9\ sinh |0| by f(\9\) and study this 
function. Its derivative f'{\9\) = 4 sin ^ sinh ^ cos -y- cosh -^(tan ^ — tanh ^-) turns into if 

|0| = 27rn, n = 0,1,2,... or 
(5.20) tan ^- = tanh 

The latter equation has infinitely many solutions \9\t, k = 0,1,2,..., where |0|o = and all 
the other roots \9\ k G {(2k - l)ir, (2k + l)ir) and \9\ k ^ 2kir, k = 1, 2, . . .. Since 

.. . . o fc o ^ \k A 9 ^ 2 \9\h \9\k \9\k 

f(0t-) = 2cosh 2cos 4cos cosh tan tanh 

■> Vl lKJ 2 2 2 2 2 2 

2 l^lfc 2 l^lfc 2 l^lfc 2 l^lfc 2 l^lfc 2 l^lfe 

= 2 cosh 2 cos 2 sin cosh 2 sinh cos = 

2 2 2 2 2 2 

and f(2irn) = 2sinh 2 (-7rn), n = 1,2, . . ., we conclude that function / has minimums at \9\ k , 
maximums at points 27rn, and it is no-negative for |0| > 0, see Figure [TJ 

The denominator <?(|0|) = sin 2 y (sinh |0| — |0|) — sinh 2 ^-(sin |0| — |0|) is an increasing func- 
tion, passing through only at |0| = 0, since the derivative g'(\9\) = ^(sinh|0| - sin|0|) + 
4 sinh 2 y sin 2 y is always positive and equals to only at |0| = 0. Hence, g is positive for 
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Figure 1. The graph of function f{9) 



Since g ^ for \9\ ^ and since f{\9\k) = 0, k = 1, 2, . . . the function /i = £ equals to only 
at points \6\k, A; = 1,2,... and is positive at all the other points \9\ / \9\k and \9\ ^ 0. The 
point \9\o = needs additional consideration. Using Taylor decomposition of function fi near 
|0|o we get that /z(|0|) = ±|0| + o(\8\ 2 ), which goes to if |d| -> 0. Therefore, /x(0) = 0. Let 

us consider the derivative = - ^^^Iqq^ 8 ^ 9 • We already know that are the 

solutions of = 0, since f(\0\k) = f'(\9\k) = 0, therefore, they are local extremal points 

of function fj,. Since fi(\9\) is nonnegative for \9\ > and /u(|#|fc) = ^j^' 1 = we conclude that 
the points \6\k are minimums of function /i. 

Moreover, since (jl is a smooth function on (0,oo), on each interval = 
0, 1, 2, . . ., it reaches its maximums, which we denote by m^. 

Let us make one more observation. Since g is a monotonically increasing function, on each 
interval \2irk, 2-7r(/c + l)], k = 1, 2, ... it has its local minimums at points 2irk. Then, taking into 
account that / has local maximums at points 2irn, n = 1, 2, . . ., we can estimate the function 
a on each of these intervals from above by the value •f( 2 ^( fc + 1 )) = 2smh (7r(fc+i)) ^ _ ± o .. 

r J g(2iTk) 7rfc-2sinh^(7rfc) ' ' ' 

From here we obtain that ^ 2smh (7r(fc+i)) _^ q w -^.j 1 ^ _^ ^ Therefore, the function 

TTfc'^t Slllrl [7VK j 

/a(\9\) decays to when \9\ goes to infinity. Note, that we can estimate tuq from above by 

/(2tt) _ 4sinh 2 7r 
3(77) 2 sinh 7r+7r cosh 7T— 37T ' 

The graph of function fi is on the Figure [2j □ 



This lemma allows us to prove the following theorem. 



Theorem 5. Let A = {x , z 1 ) be a point such that Ix 1 ^ = 0, x\ (1) < Xq(1) and uiq is a global 



maximum of fi. Suppose that 



zf(l)+ z f(l) 



-2*1 (1)^(1) 



< mo and that \9\ is a solution of the equation 



(5.21) 



zl(l)+zl{l) 



-2^(1)^(1) 



mi 



Then there exist more than one spacelike geodesies joining the origin with A. 




The equations of geodesies are given by 

(5.22, 

sinh ^ V z 

2/ n 2/ x + X l{ 1 ) ■ 2 /^l* 



and 



sm 2 f V 2 



Zl (t) = -*l±^(sinh(|^|i) - \0\t) + gi(1)+ 2 !|. (1) (sin(|6>|t) - \0\t), 
4 sinlr ^ 4 shr l -f 

(5 ' 23) *2(t) + 4(f) = X ^W^§W\ ( cosh(|0|t) - co6(|0|t) - sinh(|0|t) sin(|0|t)) 

8 sin Y Sinn Y 

x((xg(l) + x?(l))coBh(|e|t- |0|) + 2x (l)x 1 (l) sinh(|e|t - \0\) - 2sg(l) + 2x?(l)), 
-^f o n -i 3 a c A > m 0) i/ien i/iere are no geodesies of any causal type joining and A. 

Proof. Given the values (ar, z 1 ) of the end point we find \9\ as a solution to the equation (|5.2ip . 
For any of these solutions we get for t = 1 from the expression W~ 1 (i)x{t) = x(0) 

. |0| /xoOO+^iO) , -x (l)+xi(l) 

- , m M f x (l) + Xl (l) -x (l) + x 1 (l) ^ 

(5 - 24) Il(0) = T l, - 3*Ul e-M - 1 j ' 

. |0| / x 2 (l)-ix 3 (l) ^ x 2 (l)+zx 3 (l) ^ 



i 3 (0) 



i{0| ( x 2 (l) - ix 3 (l) _ x 2 {l) + ix 3 {l) \ 
2 V e~ m ~ 1 e^l - 1 )' 
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From here we calculate 

(5.25) v l = -±l + x\ 



■2 , -2 

^2 = X 2 + ^3 



X 



i) +*?(!)) 



| 2 (x 2 (l)+x 2 (l)) 



4 sinh 2 f- 



4sin 2 f 



l(isinh 2 f 



(5.26) x 2 + x 2 = LI ( e 1^1 (^o(l) + xx(l)) 2 + el e l(-x (l) + x x (l)) 2 



(5.27) x xi = ^1 (e-l e l(x (l) + xx(l)) 2 - el e l(-x (l) + xx(l)) 2 ) . 

16 sinh 2 |1 v ' 

We change in general equations the values of Sj- and by (|5.25p and get (|5.22p . 

Next, for each solution to (|5.2ip we find an orthogonal transformation (|4.9p of z-space that 
fixes zi-coordinate and therefore leave the expression z\(t) + z\{t) = z\{€) + z 2 (t) invariant. 
Therefore using (|4.1ip . we substitute there the necessary combinations by (|5.25p and (|5.26p 
and get (pT23|) . 

We observe that according to Lemma [2] projections of geodesies on (x2, X3)-plane are circles 
with center at (~ X feP > ]§\) an< ^ ra dius ^j^. The projections on (xo, xi)-plane are suitable 
brunches of hyperbolas passing through (0,0), see Lemma [TJ The parameters of the circles 
and hyperbolas can be rewritten in terms of the final point of a geodesic by making use of (|5.ip . 

Notice that if (1) + z|(l) — > 0, the quantity of solutions to (|5.21|) is growing and the values 
of solutions tends to 2irk. In this case (|5.17|) shows that the part of the velocity v% tends to 
and we come into particular situation of the Theorem [2j □ 

Remark 2. Theorem gives an answer about existence of spacelike geodesies and allows us 
estimate their cardinality. Unfortunately the estimation is not complete since we use only the 
orthogonal transformations in z-space leaving invariant zi-coordinate. The difficulty is pure 
technical, since it is complicate to find the relation between the coordinates of finite point 
and the values of \6\. In the case of orthogonal transformations in z-space leaving invariant 
zi-coordinate the relation is expressed by equation (|5.2ip . 

Remark 3. Given a point A = (x 1 ,z ) on the surface \x l \\ = 0, z 1 ^ 0, there are no timelike 
geodesies joining to A. 

5.2. Connectivity between (0, z 1 ) and (x , z 1 ), where z 1 = const. Let us use the hyper- 
spherical coordinates (r, ip, £i, £2): 

xq + ix\ = re*^ 1 cos ^, X2 + 1x3 = re l & sin ^ 



Then the horizontality conditions (|2.8p become 



Zi = --r (Cicos - + £ 2 sm -), 
z 2 = ~r 2 (<£sin(£i + £ 2 ) + sin</2COs(£i + £ 2 )(x2 - xi)), 



h = \r 2 : ((f cos(£i + £ 2 ) + sin<£sin(£i +6)(^i -^2))- 



Then, 



•2 , -2 1 -2 1 K^P 1 1 /2 2 P , ^2 ■ 2 i P\ 

z 1 +z 2 + z 3 = -r [— + Ci cos -+( 2 sin -). 
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Theorem 6. A smooth curve c{t) is horizontal with constant z- coordinates if and only if 
c(t) is a straight line in a ^-dimensional affine subspace: c(t) = (dot, . . . ,a$t, z\, z\, z\) with 
ctQ, . . . Q3 £ K 4 and a 2 , + . . . a 2 / 0. 



Proof. Let c(t) be a horizontal curve with constant vertical components. Then the equation 

h 
i 1 



lil 2 = ir 2 ^ 



+ Ci cos 2 § + H sin 2 f ) = implies 



<p = 0, £1 cos - = 0, £2 sin - = 0. 

iFvom the first equation we conclude that ip = p° is constant. In the case when ip° 7^ ir + 27r/c 
and </?° 7^ 27rn, /c, n £ Z, we see that £1 and £2 are constants £j and £2- Therefore, horizontal 
components of the curve c{t) are of the form 

n If n 09° n 03° n 09° 

xo = t cos cos — , x\ = t cos £2 sin — , X2 = t sin cos — , X3 = t sin £2 sin — . 
If if = vr + 2vr/c, ifceZ, then 

xq = 0, xi = ±tcos£2) %2 = 0, £3 = itsin^i 
and if </r = 2-7rn, n G Z, then 

xo = ±icos£i, xi = 0, X2 = ±tsin£°, X3 = 0. 

Conversely, let us assume that c(t) = (a^t, . . . , a%t, z \ , z\ , z\ ) , where 2J, . . . ,z\ are some con- 
stants. Set at = (aot, . . . , a^t). Observe, that (L\a, a) = (ha, a) = (ha, a) = for any vector 
a = (ao, . . . , 03). Then, 

li = = -(li(ai), (at)) = ^(Iia,a), z 2 = = -(I 2 (at), (at)) = ~(I 2 a,a), 

1 ■ i 

z 3 = = -(7 3 (at), (at)) = -(/3a, a). 

This implies that z-coordinates are constants. □ 

6. Reachable sets by geodeics on Lorentzian Heisenberg group 

In this section we would like to compare the results obtained for H-type Quaternion group 
with Lorentzian metric and 3-dimensional Lorentzian Heisenberg group H^. 

We remind that is a triple (R. 3 ,H, g), where M 3 is equipped with noncommutative mul- 
tiplication law 

(x ,xi,z) o (x'^x'^z 1 ) = (x + x' ,xi + x[,z + z' + ^(xix' - x x[)), 

the subbundle H is a span of two left invariant vector fields Xq = + \x\j^, X\ = — 

l x °~§z> ^ or wn ich [-^Oj-^i] = % = Jjj an d g is a Lorentzian metric on H defined by 

g(X Q ,X Q ) = -l, g(X 1 ,X 1 ) = l, g(X ,X 1 ) = 0. 

In [H] authors investigated the connectivity by geodesies in H^. In particular, they obtained 
the following result. 

Theorem 7. Let A = (x,y,z) be a point such that x > (x < 0), —x 2 + y 2 < 0, ^ B < 1. 
T/zen f/iere is a unique future directed (past directed) geodesic, joining O = (0, 0, 0) with the 
point A. Let 9 be a solution of the equation 

Az \9\/2 



-x 2 + y 2 sinh 2 (|0|/2) 



coth(|0|/2). 
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Then the equations of timelike future directed geodesic 7 : [0, 1] — > are 

\9\ \9\ \9\ \0\ \6\ 

x(t) = sinh 2 (Ut)( a; ( cot h(Ut) coth(U) - 1) + y(coth(Ui) - coth(U))), 



2 / / »v l 2 ' v 2 



y(t) = sinh 2 (^t)(y(coth(l|t) coth(^) - 1) + x(coth(^t) - coth(^))), 

|fl|t- S inh(|fl|t) 
l) " |0| - sinh(|0|) ' 

Moreover, the authors obtained the following result about the reachability by causal Hamil- 
tonian geodesies. 

Theorem 8. Let us define the following sets 

Rt = {-x 2 + y 2 < 0, -i!£L<l}, 
x z — y A 



R sp = {-x 2 + y 2 >0, _ 2 " M ,. 9l < 1}, 



4\z\ 
-x 2 + y 2 l 
Ri = {-x 2 + y 2 = 0, z = 0}. 

Then there exists a unique geodesic connecting the point O = (0, 0, 0) with a point P that 
belongs to one of the sets Rt, R sp or Ri. Particularly, if P £ R t , then the geodesic is timelike, 
if P € R sp , then the geodesic is spaselike, and if P £ Ri, then the geodesic is lightlike. 

The connectivity on is dependent on the solutions \9\ of the equation 

4z 



-x z + y z 



where the function /i(t) = sinl ^y 2 ) — coth(r/2) is strictly decreasing on the interval (—00,00) 

' > 1, i. e. if 
le origin with 



4z 



-x 2 -\-y 2 



from —1 to 1. It means that if the point A = (x,y,z) is such that 

A £ (R t U R sp U Ri) c , then there are no geodesies of any causal type joining t 
A. In particular, there are no Hamiltonian geodesies joining the origin with the points of the 
surfaces {7^1 = 1, -x 2 + y 2 < 1} and { .Jj^ = 1, -x 2 + y 2 > 0}. We observe that in [6] it 
is shown that it is possible to find nonHamiltonian lightlike geodesies joining with the points 
of these surfaces. 

Thus, in contrast with Ql, the Heisenberg Lorentzian group has the property of uniqueness 
of geodesies starting from the origin with given tangent vector. It happens due to lower 
dimension of the "spacelike" part of H^. 
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